L 00 (L 00 )-BOUNDEDNESS OF DG(p)-SOLUTIONS FOR NONLINEAR 
CONSERVATION LAWS WITH BOUNDARY CONDITIONS 

LUTZ ANGERMANN AND CHRISTIAN HENKE 

ABSTRACT. We prove the L°° (L°° )-boundedness of a higher-order shock-cap- 
turing streamline-diffusion DG-method based on polynomials of degree p > 
for general scalar conservation laws. The estimate is given for the case of several 
space dimensions and for conservation laws with initial and boundary conditions. 



1. Introduction 

In this paper we extend the analysis of a shock-capturing streamline-diffusion 
DG-method for hyperbolic conservation laws in several space dimensions which 
goes back to [JJS95]. The original DG-method from [JJS95] is based on polyno- 
mials of maximal degree p > (DG(p)-method) and is applied to a pure Cauchy 
problem. Here we formulate the method for scalar conservation laws with initial 
and boundary conditions. To describe the further features of the method, we recall 
the following sufficient conditions for convergence of a sequence of approximate 
solutions [Sze89b, Remark 1.2]: 

(1) uniform boundedness in the L 00 (L 00 )-norm, i.e. L°° in time and L°° in 
space, 

(2) weak consistency with all entropy inequalities, 

(3) strong consistency with the initial condition. 

In the case of an unbounded domain, the condition (1) can be replaced by the 

(1*) uniform boundedness in the L°°(L 2 )-norm, 

which was done in [JJS95]. To the best of our knowledge, the only attempt to 
prove the L 00 ^ 00 ) -boundedness without using a finer auxiliary triangulation is 
given in [JSH90] for the case p = 1. This proof can be extended for p > 1 if the 
shock-capturing terms are defined on finer triangulations [Sze91]. Thanks to (1*) 
this is not necessary for the DG(j?)-mefhod in [JJS95]. Our result presented here 
uses the skeleton of the proof from [Sze91], which is based on choosing the test 
functions v = I^U^ 1 ) with a large even number q, where 1^ is the Lagrange 
interpolation operator and U denotes the approximate solution. Within this proof 
we use a new algebraic argument to verify the coercivity of the shock-capturing 
term when v = I^(U q ~ 1 ). 

Let us recall the key points of the DG(p)-method. First, we have to choose 
the numerical flux on the element boundaries. In contrast to [JJS95], where a 
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strictly monotone numerical flux is necessary, we may also use monotone nu- 
merical fluxes such as Engquist-Osher fluxes. Second, there are two stabiliza- 
tion mechanisms. The DG(j?)-method from the last mentioned reference contains 
a streamline-diffusion term and a residual-based shock-capturing term. The first 
term adds an anisotropic artificial viscosity and the second one introduces some 
isotropic artificial viscosity locally where the solution is nonsmooth. 

The paper is organized as follows: In Section 2, we prepare some basic material 
on scalar hyperbolic conservation laws with initial and boundary conditions. Then, 
in Section 3 we introduce the DG(p)-method under consideration, and in Section 4 
the condition (1*) is verified. After this we present our main Theorem 5.1 which 
is proved in Sections 5 and 6. Here we give some background material on spatial 
and algebraic numerical ranges and extend the condition (1) to the case p > 1. 

The entropy consistency and the consistency with the initial condition will be 
proved in a forthcoming paper. 

2. Hyperbolic conservation laws with Boundary Conditions 

Let Qt = (0, T) x CI C R d+1 , T > 0, d € N, be an open time-space domain 
with boundary T>t = (0, T) x T, T = dCl and with outward unit normal n. In this 
time-space domain a point with position x = (x±,X2, • • • , x<i) T at time t = xq has 
the coordinates X = (xq, x) T . Standard notation is used for the space of functions 
of bounded variations BV(Qt), Lebesgues spaces L q (Qt) and Sobolev spaces 
W l > q (Q T ), I € N, 1 < q < oo. 

We consider for u : Qt — > K the initial-boundary value problem 



with the following boundary condition: For all k € R, r € 

(sign(7u(r) - k) - sign(g D (r) - k)) {f{-yu(r)) - f(k)) ■ n(r) > 0, (2.3) 

where T = (•, f) T : R -> R d+1 , u : O -> R, g D : S T -> R are given smooth 
functions and 7 : Qt — > Xt denotes a trace operator. The function sign : R — > R 
is defined by 



Due to the hyperbolic nature of (2.1), a boundary condition of the form u = gr> on 
St usually over-determines the problem. The generalization of the inflow bound- 
ary condition (where f'(gD) - n < 0) for nonlinear / also leads to a problem that is 
not well-posed. This difficulty does not occur in (2.1) — (2.3), because the solution 

u e of 



L(u) = V • T(u) = Oin Q T 
u(0, •) = uq on Q 



(2.1) 
(2.2) 




x = 0. 



eAu e + L(u £ ) = in Qt, 
u e = g tD on S T 
u e (0, •) = u t o onfl, 



(2.4) 
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converges a.e. to a function u G BV(Qt) as e — > , which satisfies (2.1) — (2.3) 
[B1RN79, Theorem 1]. Moreover, we can use this vanishing- viscosity method even 
in the space L°°(Qt)- It is possible to define a well-posed initial-boundary value 
problem, which admits a unique weak entropy solution u G L°°(Qt) [MNRR96, 
Definition 7.2, Theorems 7.28, 8.20]. 

Let us shortly recall the concept of entropy pairs. We say that Q = (77, qi , . . . , q^) 
is an entropy pair if 77 : E — > E is continuous and convex, the entropy flux 
q 3 : : E — > E is continuous and 77, gi, . . . , qj, satisfy for all u 6 E the compati- 
bility condition 

V '(u)f' ] (u)=q' J (u). (2.5) 

For scalar conservation laws this is trivially fulfilled if the entropy flux is defined 

as 

<&(«) = / v'^f'^dr, l<j<d. (2.6) 

3. Formulation of the DG(p)-method 

In this section we introduce the DG(y?)-method. To discretize (2.1) — (2.3), let 
Qn,n+i = (t n , t n+ i) x f2, Q n = {t n } x for the sequence of discrete time levels 
= to < t\ < ■ ■ ■ < ij\r, iV G N, be a time-space decomposition of Qt- The 
boundary is defined by S nin+ i = (t n , t n+ i) x T and E n = {t n } x T. 

Consider an affine decomposition Tf^ of Q n ,n+i belonging to a family of quasi- 
uniform, admissible decompositions of Q n n +i, cf. [EG04, Definition 1.49, 1.53, 
1.140] into simplices or quadrilaterals T and write Th = Un>o • L et be the 
diameter of T and h the maximal diameter of all T G T£. 

Set 

W£ = {w G L 2 (Q n , n+1 ) : u/| T G P P (T) VT G 7^} , = J] (3-1) 

n>0 

where 

P P (T) = span {X a } ,XeT 

is the space of polynomials of maximal degree p defined on T. 

We are now ready to define the DG(p)-method and introduce by R l n n+1 , R l n the 
set of all interior faces of Q n ,n+i, Qn and by A n n+ i, A n the set of all boundary 
faces. We further set R n ,n+i = R n n+i^^-n,n+i and R n = i?^UA n . In order to be 
able to describe discontinuous functions, we denote by r the common face shared 
by the elements T + = T and T~ . We also define the normal vectors = and 
nTp on r. Then we introduce the notation 

v ± (x) = lim v(x — fm^), v±(x) = v(t n ± 0, X\, . . . , xa) (3.2) 

and 

{v} = -(v + + v~), [vn] = -u + n + + v~n~, fvj = v + — v~ . (3.3) 
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If v is a function on A n n+ i or A n , we write v~ = By C we will denote a 
positive constant independent of h, not necessary the same at each occurrence. 
Introducing the bivariate form 

a(v,w) = J2 2 {/ L (p)wdX + [ {J r {v)-F{v + )-n + )w + ds 

n=0 T<=T£ ^ T 9T 

(3.4) 

the DG(p)-method for (2.1) — (2.3) can be formulated now: Find U € Wh such 
that for n = 0, 1, • • • ,N,U= U Qn n+1 e W% and for all v £ Wff 

a(U,v)+ {{6L(U),T'(U)-Vv) 0T + e(VU,Vv) OT } = 0, (3.5) 
TeT% 

where 

5 = ( 5(C/) = C 1 / l T(||^ / (f/)||p)" 1 , 

e = i(U) = max f C 2 h 2 ~P R(U), C 3 h p+l / 2 ) , < < -, 



12(C0|t = ma X (|L(C/)|) + -L (max (|[^(l7)n]|) + max (C T M) 

d*T = {x E dT : x i R n+1 ) 

and C%, C2, C3 > 0. As mentioned in the introduction, (3.5) contains a streamline- 
diffusion term and a residual -based shock-capturing term. Due to the /i-dependency 
of this term, the DG-method can be interpreted as a discrete vanishing-viscosity 
method. 

The numerical flux JF{U) is given by 

P{U) = {F(U)} ■ n+ + C T (U+, U-,n+)lV} (3.6) 

and 

i n+ = ±(l,0,...,0), 

C T (v + ,v-,n+) = { C 9Q > £ \T'{v~ + S [vJ) ■ n+\ ds v~ = g D , 



C ° > i Jq \T'(v~ + sfuj) ■ n+\ ds otherwise. 



(3.7) 

Within this framework there are the following well-known numerical fluxes: 
The Engquist-Osher flux if /(0) = 0: 

Cri^-, „♦) . \ C X -$?< r + *« ■ " +l * \ = (3.8) 
I Cq = ^Cq otherwise. 

The Lax-Friedrichs flux: 

= sup \F'{z)-n + \ v~=g D , 
C T (v + ,v-,n + ) = < ze[v+,v-] (3.9) 

[C$ = ±C$ Q otherwise. 
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Further Co is a positive constant satisfying 

C > ||^'||o,oo,r = max\\P{x)\\p (3.10) 

and 

C™,C$<C . (3.11) 

4. Preliminaries 

In this section we want to verify the uniform L°°(L 2 ) -bound (1*) which is based 
on [JJS95]. As a byproduct of the study of this technical result, we introduce 
some notation and prepare the basic material for the L 00 (L 00 )-bound presented in 
Section 5. To make this precise, we choose T){U) = U 2 /2, ip = 1 and v = r)'(U)ip. 
The main theorem of this paper can be obtained by applying the test function v = 
/£(?/(£/», where r){U) = U q /q, (p = 1 and l v h : C(Qt) -> W h is the Lagrange 
interpolation operator. Here, due to the fact that rj'(U) ^ Wh, an interpolation 
or projection operator is necessary. Notice that in this case we get an additional 
difficulty to estimate terms which contain the difference rf{U)ip — l?(?/(f7)<£>). 

By the definition of the bivariate form 

N-1 

b(v,w)=a(v,w) + Y J E {S(v)L(v),F'(v)-Vw) QT (4.1) 

n=0 T€T,l l 

we have, for w = rj'(U)(p and an entropy pair (7?, q) satisfying (2.5) and (2.6), that 

N-1 r 

b(v, v '(v)Lp) = J2 E (W«)^(iV(#) P )) 0T 



n=0 TeT h " 



+ / &(v)nW(v)<pds+ [ QrMr/(v)<pds 
JdT z JdT 

+ / Q(v)-mpds- I Q(y)-V(pdX. 

JdT J\0.tM]xQ 



(4.2) 



I&T J[0,t N ]xCl 

As usual in DG-methods we consider the different behaviour of inner and boundary 
faces 

b(v,rf(v)<p)= [ i](v^)^ N dx- [ rj{v Q _ )<p° dx 

-[ Q(v)-V V dX 
J[o,t N ]xn 

5 



i=0 

where 

N-1 



E (f,r),v,<p) = J2 E {5(v)L(v),r(v)-V( V '(v) V )) QT , (4.4) 
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E 1 (f, n, v, <p) = ]T / (r,(yl) - vK) - r/(vl){vl - <)) V n dx, (4.5) 

, n JO, 



N-l 



n=0 
N-l 



E 2 {f,r 1 ,v,y) = Y J E / {lQ(vM-lHvM{v'(v)})^ds, (4.6) 



N-l 



E 3 (f, v ,v,<p) = J2 E fc$[Tf'(v- + rM)drM 2 <pd8, (4.7) 



£ 4 (/,r?,^) = E E {»/(«)}) (4-8) 

n=0 T6A n , n+ i T 

AT-1 ! 

E 5 (f, V ,v,<p) = Y, E / C 'o n / ^"fe + rH)rfrHVd S , (4.9) 

„-n ^^A . . ■'r JO 



n=0 rGA„, n+ 
AT-1 



n=0 rGA„,„ +1 • 7t V ' 

(4.10) 

Next, we will show the nonnegativity of Ei (f,r), v, ip). By the convexity of rj, 
this is true for E\ . In order to treat E^ , we consider the expression 



(q(v+) - G(rT) - ^ (^ + ) - .F(tT)) (r/(*, + ) + rf{v~)) 
(jf {Q' -J r 'v')dr + £ T' (r{ - ^ {r]' (v+) + r,' (v'))^ drj 



n 



T\v~ + s\v\) • n+ ( 7]'(v- + sH) - - (r/(i; + ) + i/(t;-)) [u] ds. 



(2.5) Jo 

By the properties of convex functions, it follows that 

(»/(«- + «H)-»/V)) H >o 

and 

(t/(«- + sH) - r/V)) H = (,/(«+ - (1 - S )H) - H < 0. 

Then we have 

v '(v- + slv})-U V '(v + ) + V '(v~))) [v 



\ { V '(v- + sM) - V '(v-)) + \ { V '(v- + aM) - rf(v+)) ) [v\ 



< + «H) - »?V)) HI + o IOtV + S H) - »/V)) H| 
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\ ((»/(«" + a[v}) - T,'(v-j) - {v'(v- + slv]) - rf(v + ))) [«] 



1 , „ ^ _^ , , 1 



j/(« + )-7/(0)N = o / V"(v-+rlv])drlv] 2 , 



2 

which immediately implies that 



o 



(q(v + ) - Q(«-) - \ (T(v + ) ~ JP(V~)) (»/(«+) + 7?' («")) 

< / IJ 7 '^" + -n+ 



o 



r,'(v- + 8[v})-Ur/(v + ) + r/(v-)))lvl 



ds 



1 



^/ l^'iv" + slvj) -n + \ I rf'{v- + r\v\) dr\vf ds 



o 



^\ I \F'{v~ + s{v\)-n + \ds I' r]"{v- +r{v\)dr{vf. 

Finally, having in mind the fact that rf' > 0, we can use the definitions (4.7) and 
(3.7) 

E2(f,V,v,ip)+E 3 (f,r],v,(p) 

* E ^(^-^V / (^+^H)^ + |^)(^\ // (^+rH)^)HV^ 



n=0 



n,n-\-l 

>0, 

(4.11) 

where we have used nonnegative test functions ip. The same arguments as before 
lead to E 4 (f, rj, v, if) + E 5 (f, rj, v, tp) > 0. 

Remark 4.1. The local condition (3.7) 

C%-\] \T , {v-+$M)-n+\ d£>0 

for (4.11) allows a smaller constant Cq than [JJS95, Remark 2.5] 

n n l,,-,,, Jo I[ (n"(v- + rjv}) + r/>+ - r M)) dr ds 

C ~ olK IIo,oo,r 73— — — > U. 

Jo V \ v + r i v l) dr 

Thus, the corresponding numerical flux J r (U) is a monotone flux function (f7 + h4 
•F(£/) is increasing and J7~ h-» -F([/) is decreasing). As mentioned in [JJS95], the 
requirement that F(U) is a strictly monotone numerical flux, e.g. 

C o ~ \ I l^iv- + £M) • n+| ^ > e > 0, 
2 JO 
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is necessary for (2). More precisely, the condition 

JV-l 



E E / hR{U) 2 dX<C\\u \\l u 



n=0 T&Ti 



has to be fulfilled. 

Let us now consider the equation (3.5) with v = r)'(U)tp: 



N-l 

b(U, V '(U)<p) + J2 E ^)(W,V(r/(C%)) 0j , 

n=0 T£T/™ 

JV-l 

= / V (U?)<p"dx- / ^(C/^^dx+E E e(£/)(W,V(r/(£/>)) 0)T 

J!) n=Q TgT^ 1 

5 

- / Q(U) -VipdX + Y] Ei(f, 77, 17, ¥>) - F(/, r?, U, = 0. 

J[o,t N ]xn i=Q 

(4.12) 

Therefore, in the case rf(U) = U 2 /2 and = 1, we obtain 

JV-l 



\ f (U») 2 dx+J2 E W(W,W) 0jT 

5 

+ E U 2 /2, U, 1) = i /" (u ) 2 cfe + F(/, f/ 2 /2, 17, 1), 



i=0 

where we have used = Uq. By the help of Young's inequality we see that 

JV-l 



F(f,U 2 /2,U,l)<Y. E J T \c^\Uf + lc^g 2 D ds 



n=0 t6A„,„ + i 



(4.14) 



U'/2,U,l)+F 2 (f,gf /2,g D ,l), 



thus we arrive at F 4 (/, U 2 /2, U, 1) + E 5 (f, U 2 /2, U, 1) - Fi(/, U 2 /2, U, 1) > 
and 



JV-l 



^IllUn + E E H*(t0 1/a ^(EOIIo A T < |ll^o||g j2 ,o + j^olMloW 

n=0 TaT£ 

(4.15) 

Finally, using for £at_i < i < ijv and g £ N the identity 

r^V 1 1 <j „ / / TTq-lr + i 



ll^(*.-)llg, ff , n = ll^llg igi n-« J J^-\t',x)dwT(U(t',x))dxdt', 



(4.16) 
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Young's inequality yields 



l|E/(Mfen < IP-*||g, 2 , n + 2Ci ^ " IP«',)ll§ A n*' 

e/ft (4.17) 

+ 2 £ ||*(E0 1/a £(tOII§ A T 
r67 ^-i 

and a Gronwall argument estimates the right-hand side of (4.17) by means of the 
left-hand side of (4.15). The quasi-uniformity of jv^" 1 j ensures the bound- 
edness of the Gronwall constant and the following theorem results. 

Theorem 4.2. Let Q be a domain with a Lipschitz boundary and {Th}h>o be a 
quasi-uniform family of decompositions 0/(0, T) x O. Let U be a solution of (3.5) 
satisfying the assumptions (3.10) and (3.11). Then there exists a constant C > 
independent of h, such that, for all t € (ijv-lj£iv)j 

\\U(t, .)|| ,2,n < C(||«o||o,2,n + ||&d||o,2,s t ) • (4.18) 



5. L 00 (L 00 )-BOUNDEDNESS OF DG(p)-SOLUTIONS 

As announced above, in this section we prove that ||£7||o,oo,Q T is uniformly 
bounded. The main idea is to control the interpolation error U q ~ l — /^(L^ 9-1 ) in 
the second argument of the bivariate form by the aid of the special shock-capturing 
term. At the end of this section we formulate two corollaries which are conse- 
quences of the limiting process h — > and of the special case p = 0, respectively. 
Our main result is the following. 

Theorem 5.1. Let Q be a domain with a Lipschitz boundary and {7h} ft> g be a 
quasi-uniform family of decompositions of(0, T) x Q. Let U be a solution of (3.5) 
satisfying the assumptions (3.10) and (3.11). Then there exists a constant C > 
independent of h such that 

\\U\\o,oo,Q T < C (||uo||o,oo,n + ||sd||o,oo,£ t + 1) • (5-1) 



The proof is based on the next lemma which contains the extension of [Sze91, 
Lemma 3.3] and [Sze89a, Lemma 4.2] for p > 1 and which is proved in the last 
section. 

Lemma 5.2. For Lagrange finite elements with a shape regular family of meshes 
{Tfr} h>0 there is a constant C > independent ofq and h such that for all v £ Wh 
and q = 2m, m G N : 

{Vv,Vll(v^ l )) 0T >C f \\Vv\\ 2 P \\v\\$£ T dX, VT€T h n . (5.2) 

Jt 
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Proof of Theorem 5.1. Setting v = I?(r)'(U)(p) in (3.5) with 77(f) = v q /q, <p = 1 
and q > 2 an even natural number, we obtain 

5 



- / {u^y dx + Y,E, t (f,uy q ,u,i) - {biu,^- 1 ) -biujliu^- 1 ))) 

q J ° i=0 
N-l 

+ E E e(U){VU,VI p h (Ui- 1 )) QT = - / (u o ydx + F(f,Uy q ,U,l). 

n=0 TeT^ 1 n 



(5.3) 

Here the key point is that the interpolation error is bounded by the isotropic shock- 
capturing term. To see this, a careful consideration of the interpolation operator 1^ 
is necessary. 
By 

/ hj c (U)n} + C T lU}dx = 0, 

J Rn+l 

we express the interpolation error as 

(b^u^-KujKu^ 1 ))) 

N-l r 

r h E E {/ V-HU)(U^-Il(U^)) dX 



(4-1) n=0 TeT% 



+ j T 5{U)L{U)F(U) ■ (VU*- 1 - VII (U*- 1 )) dX (5 4) 

Jd*T z 

r N-l 4 

+ / C T \U\ (U^-II{U^)) ds} = E E^t- 

Since U\t S Pp(T'), we deduce that |C/| p +i,oo,t = 0. Arguing as in [Sze91][p. 
765], we may write for q > 3 

|tf*-Vi,oo,T < ^ P+1 /i~ P+1 ||V[/||g i0O , T ||C/||^ iT . (5.5) 

Together with a standard interpolation error, we thus conclude from the first term 
of (5.4) that 

K| < W-iIW^w^t -T(u)\dx 

< C^ +1 |^- 1 | p+1 , 00 , r J \V -T(U)\dX 

< c q v +l hl\\vu\\l^ T \\u\\lt T I \V-T{U)\dX 

(5.5) Jt 

< Cq? +1 h 2 T [ \\VU\\l^ T \\U\\l* T \V-T(U)\dX 

JTn{\u\>i} 
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+ Cq p+1 h 2 T [ ||V?7||g i00|T [V • F(U)\dX 

JTC\{\U\<1} 

< CqP +1 h 2 T mM\V-T(U)\)\\U\\ g ^ T [ ||Vl7||iU r d# 

T JT 

+ Cq p+1 h 2 T ma^{\V ■ F{U)\) / || VU\\l ^ T dX. 

T JT 

By the quasi-uniformity of {T^} h>0 and an inverse inequality, we obtain that 

[ ||Vi/||jj T dX < C [ \\Vv\$dX. (5.6) 

JT ' JT 

Thus we have, by Lemma 5.2, 

\A l T \ < C/ +1 4max(|V ■ T{U)\) { (W, V ^(U^ 1 )) Q>T + ||VZ7||g i2 , T } . 

In a similar fashion, we can estimate the complete right-hand side of (5.4). Conse- 
quently, by (4.15) we conclude that 

\{b(U,U^ l )-b{UJl(U^)))\ 

N-l 

< Cq p+1 Y Y h TR(U) (W, VlfXUi- 1 ))^ + Ch^ +1 . (5.7) 



n=0 TeT h 
Inserting this into (5.3) we obtain 

JV-l 



0,T 



/ (u N ) q dx - c q p+2 Y Y h rR(u) (w, vif l (u^)) T 

Jn n=0 Te% 

N-l 

+ E KU){VU,Vll{U^ 1 )) ()T +{5(U)L(U),F'{U)V(U^ 1 )) 

5 

+ Y Ei(J, U q /q, U,l)< / (u ) q dx + qF(f, U q /q, U, 1) + Ch P T q p+2 . 

(5.8) 

To proceed with the treatment of interpolation error, it is necessary to require that 
Cq p+2 < h" 13 , where < < 1/2. The upper limit of (3 is introduced due to 
convergence reasons, cf. [JJS95, Lemma 3.2]. However, this restriction on q does 
not prevent us to finish this proof by letting q — > oo. 
Moreover, we have 

. N-l 

/ (U N ) q dx + qY Y (W)L{U)^'{U)V{Ut- l )) QT 
5 

+ Y E i(f,U q /q,U,l)< / (u ) q dx + qF(f,Uyq,U,l)+C. 
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Note that 



E 4 (f, uy g ,u,i) + E 5 (f,uy q ,u,i) 

N-l , , , fl 



* E E I ( C o U -l [ \F'(9D + f [CI) • n + | df) [t^lPI cfa, 

„-n , c » ..•'t\ z JO / 



n=0 reA„, n+ i 

and thus, by Young's inequality, 

F(f, uy q , u, i) 



N-l 

< 

n=0 reA n , n+ i 
JV-1 



E E I (\ l\^(9D + aU\)-n + \di + C d A\lUU D l \ds 

„-n , Jr \ z JO J 



^ E E \\c^W X W\** + ^(3-1^3^2^^ 
(3.10) ^fc^ 2 2 

= F 1 (f,Uyq,U ) l)+F 2 (f ) g q D /q > g D) l)- (5.10) 

So we obtain 

£ 4 (/, C/V9, U, 1) + £ 5 (/, 17V«Z, U, 1) - Fi(/, 179/g, U, 1) > 0. 
Altogether we get that 

„ N-l 

/ (U») q dx + qJ2 E (^(CO.^CCOVCL^ 1 ))^ 
n=o re7^ 

< / ( no ) 9 (irE + Co( (? -l)(2g)-^3^|| to ||^ ET + a 

(5.11) 

By repeating the arguments given at the end of the previous section, we summarize 
that 

sup||Ef(V)||o,,,n < C\ f|K|| 0>g ,n + 3i^ [c m )\ |M| , ? ,s T + c\\ (5.12) 

g 

for 4 < g < C7i p+ 2 . Finally, using an inverse inequality we have that 



||f/|| ,oo,Q T < (Cg/i- 1 )' sup||Z7(t,-)||o, gj Q < (c/ J - 1 -^)'sup||C/(t,-)||o 



t>0 7 t>0 



= C* exp [Cdcf x In (/t" 1 )) sup||C/(t, OHo^n- 

t>o 

Setting q = Ch p+ 2 we get 

||C/||o,oo,Q T < C d ^ /(P+2) exp (C^/(P+ 2 ) In (/i" 1 )) sup||C/(i, .)||o,,,f 



<C7sup||l7(t,-)llo, g ,n 

i>0 



*-° (5.13) 



in the case of h < 1 , which concludes the proof. □ 
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Corollary 5.3. Under the assumptions of Theorem 5.1, the estimate 

||^||o,oo, Q T — ll u o||o,oo,n + HffD |[o,oo,E T + 1 (5-14) 

holds for h — > 0. 

Proof Using (5.12) and (5.13) with h = Cq"^ 2 ^ 13 -> 0, the statement imme- 
diately follows. □ 

Remark 5.4. The L 00 (L 00 )-boundedness of DG(p)-solutions withp > 1 was also 
considered in [Sze91]. However, since an inequality of the form (5.2) was proved 
only for the case p = 1, the shock-capturing term was realized on finer auxiliary 
triangulations using polynomials of first degree. Hence, the bound C q q p+2 < h~P 
is necessary, which is true for q < Cln(l//i), cf. [Sze91, (3.16)]. This gives 

||C/||o,oo,Q T < Csup||£/(i,-)Ho,<z,fi, 
t>o 

for h — > 0, where 1. 

Corollary 5.5. Under the assumptions of Theorem 5.1, the estimate 

||^||o,oo,q t < ||«o||o,oo,n + ||5d||o,oo,s t V/i > (5.15) 

holds for p = 0. 

Proof Let p = 0. Then we get (biU,^" 1 ) - b(U, if^U^ 1 ))) = 0. Conse- 

§_ 

quently, there is no need for the bound q < Ch p+ 2 , and we can conclude with 
letting q -)• oo in (5.12) and (5.13). □ 



6. Proof of Lemma 5.2 

Until now there is no proof of an inequality like (5.2) for q ^ 2. To the best of 
our knowledge only special cases for linear ansatz functions on triangles respec- 
tively tetrahedrons are available, cf. [Sze89a, Lemma 4.2] and [Sze91, Lemma 
3.3]. Moreover, the constant in these references depends on q. 

Using the theory of numerical ranges for bounded linear operators in Banach 
spaces, we are able to prove this inequality under rather weak assumptions. More 
precisely, the local stiffness matrix of the shock-capturing term has to be symmetric 
positively definite and an eigenvector (1, • • • , 1) T with an unique eigenvalue zero. 

First of all we need some further notation and definitions about the numerical 
range. 

Definition 6.1. Let (X, || • ||) be a normed vector space, let S(X) be the unit sphere 
and denote by X' the dual space of X. For each bounded linear operator A on X, 

W(A, \\-\\) = {f(Ax): (x,f) en} (6.1) 

with II = {(x,f) e S(X) x S(X') : f(x) = 1} is called the spatial numerical 
range. 
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Remark 6.2. 

(1) Note that the following definition is equivalent to (6.1): 

W(A, || • ||) = {f(Ax) : f(x) = \\x\\\\f\\ = 1} , (6.2) 

since 

||/||= sup |/(x)| < ||/||^=1. 
\\ x \\ 1 =1 

(2) In contrast to the spectrum o~(A), the spatial numerical range W(^4, || • ||) 
depends on the norm || • ||. 

(3) Let l qi (n) = (M n , || • Wm) be the normed vector space. Due to f(x) = 
Ya=i x if( e i) = xT Uf an( l l q2 (n) for 1/qi + l/q2 = 1, which is norm- 
isomorphic to l qi (n)', the identity ||/|| = ||y/||;?2 is valid. Therefore we 
get 

W(A, || • || m ) = {x T Ay f : x T y f = \\x\\ m \\y f \\ m = l} , (6.3) 

which is the identity case of Holder's inequality, cf. [Bau62]. In the case 
qi = 2, VF(^4, || • = W(A) is the numerical range in a Hilbert space 
due to Toeplitz [Toe 18]. 

(4) Unlike to W(A), the spatial numerical range is not necessary convex, cf. 
[NS64, S. 357]. 

If we interpret the Matrix A as an element of a normed algebra A with an identity 
element, we can define a second numerical range. For further details we refer to 
[BD71, S. 15]. 

Definition 6.3. Let A be a normed algebra, S(A) = {x £ A : ||x|| = 1} the unit 
sphere and A' the dual space of A. For x G A let 

D(A,x) = {feA':f(x) = l = \\f\\}. (6.4) 

We define the algebraic numerical range by 

V A (a, || • ||) = U {V A (a,x, || • ||) : x G S(A)} , (6.5) 

where 

V A (a, x, || • ||) = {f(ax) : / G D(A, x)} . (6.6) 

Notice that for the algebraic numerical range it is sufficient to consider only the 
identity element. 

Lemma 6.4. 

VaM ■ ||) = V A (a,l,\\ ■ II), a G A. 
Proof. [BD71, Lemma 2.2]. □ 

Next, we recall two well-known results about numerical ranges and the connec- 
tion to the spectrum a (A). 

Lemma 6.5. 

convW(AIHI) = ^(a,HI)- (6-7) 
Proof. [BD71,S. 84]. □ 
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Theorem 6.6 (Vidav). Let a £ A be a Herniitian element, i.e., V^(a, \\ ■ ||) C R. 
Then we have 

convcj(o) = V A (a, \\ ■ ||). (6.8) 
Proof. [BD71, Corollary 5.11]. □ 

Further we have a corollary which will help us to prove Lemma 5.2. 
Corollary 6.7. Let Abe a symmetric, positively semidefinite matrix. Then we have 

W(A,\\-\\)Q[0,\ max (A)}, (6.9) 
where A max (yl) denotes the largest eigenvalue of A. 

Proof of Lemma 5.2. First, consider (5.2) on the reference element. Obviously, the 
inequality is valid for v = const . Let v ^ const be given. Consider a decomposi- 
tion of Q P (T) = V°(f) © V{f) with 

V°(T) = {veQ p (f):v = const} 



{v € Q p (f ) : / VwVvdx = 0, w £ Q p (f)}. 
Jf 



Let Vj\f and V$ , Vjy- © VV = I"" be the coefficient spaces. Therefore, due to the 
definition of Lagrange finite elements we have that 

V# = span{(l, . . . , if}, dimV# = 1. (6.10) 
Let N denote the Lagrange nodes and n& f the number of degrees of freedom. 
Moreover, defining = (Vipi, . . . , V(^ ndof ) T , vjj = (v(x)) xeA f, vfc 

= (v q ~ l (x)) x& N, we have that 

W—TTq = \ £l = V M Av q N 

if Vj^-Vj^ 1 = 1. Notice that, due to the homogeneity of the quotient, such a norming 
is always possible. Using the fact that A is a symmetric, positively semidefinite 
matrix and 1 = vjj-v^ 1 = ||w IWI^a?" ]|^a/c«— i) , we obtain 

vlAvfc 1 e W(A, || • C [0,A max (A)]. 

(6.9) 

Now, it is natural to ask whether vJfAvJj 1 is bounded from zero independent of q. 
To see this, let us suppose that the eigenvalues are ordered in increasing manner 

= Ai < A 2 < • • • < A„ dof . 

Since V$ is the eigenspace of Ai, we have 

VfrLVM, (6.11) 

and we can write A in terms of a sum of dyadic products of eigenvectors £j , 1 < 

i < n d0 f 

"dof "dof "dof 

v TfAv q x l = vJf^XMjv^ 1 =J2 XiV %ftei v j7 1 =^Z X i v M^I v M 1 - 

i=\ i=l i=2 
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Again, the inclusion (6.9) yields vjf^fvj^ € [0, 1]. On the other hand, suppos- 
ing vj\f ^ const, we obtain 

vjV, v q j^ 1 -L span{^i} = V$. (6.12) 

Then, the boundedness of vJ^AvJ^ 1 , i.e. 

"dof 

implies that 

(V^VI^" 1 ))^ 

n no ~ — -^2, v ^ const . (6.13) 

\\ V M\l9 

Now, standard estimates give the proof for the reference element 

IHCt W V <2,f $ A max (i)||^- 2 || 0jOOif ||wll? 2 

ll^ll?oo 2 ||wll?2 

< A m ax(^)Apiiwiii'<r 2 iiwii^ 

< (^dof) 1 ^ 2 ^ A max (A)Ap||uAf||^ 

<n dot ^^A p (VvyF h (v^)) of , 

where A p = \\ YTi=[ I ^1 1 1 o oo f * s tne Lebesgue constant. 

Finally, we want to show the result for an affine decomposition {T^} h>Q with 

F T : f 3 x ^ J T x + b T = x € T VT € {T£} h>0 

and therefore 

Vu(x) = J^ T Vu(x), u = uoF T . (6.14) 
Due to the spectral decomposition of 



i=i 

we have that 



V N (J (y<Pj) T KV<Pi\ det( J T )\ dx 
I det ( j t)\ ^2 (Hvjf ( I (Vifj^^iipfVipi dx J vjj 



1=1 ' i,3 



> | det(J T )|^ min t; ArJ 4fL- 
(6.9) 

= |det(J T )|||J T ||- 2 2 ^i^- 1 
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and 

\Vv\* Af > |det(J T )|- 1 ||J T 1 ||- 2 |V V |2 2iT . 
Now, the shape regularity property 

II 1 1 II 1 1 U 2 — C 
completes the proof. □ 

7. Summary 

In this paper, we considered a DG-method based on polynomials of degree p > 
for hyperbolic scalar conservation laws. This method was introduced for the pure 
Cauchy problem in [JJS95]. We extended the formulation for hyperbolic conser- 
vation laws with initial and boundary conditions. Moreover, we presented a proof 
of the uniform boundedness of the discrete solution in the L 00 (L 00 )-norm. The 
analysis is based on arguments demonstrated in [Sze91] which are valid for p = 1. 
It turned out that the use of numerical ranges for bounded linear operators in Ba- 
nach spaces allows to generalize this result to the case p > 1. Future work will 
be devoted to the convergence of the DG-method for the initial-boundary value 
problem. 
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